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In the limit of infinitely heavy /b-quarks we compute the SU (3)-breaking ratio S, of neutral B- 
meson mixing matrix elements. We also present results for the ratio of decay constants /bJ/bj- 
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action with link-smearing for the fe-quarks, and the Iwasaki gauge action. Here we report on 
our results from the 16 3 x 32 x 16 ensemble (a -1 = 1.729(28) GeV) which we use to verify our 
method. We improve upon our earlier work by including 0(a s pa) matching for the computation 
of the decay constants and extrapolating to the physical point using chiral perturbation theory. 
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1. Introduction 



Neutral B-B-mixing is a phenomenologically important quantity because it allows us to de- 
termine CKM matrix elements. Within the Standard Model the dominant contribution is given by 
box diagrams with top quarks as shown in Figure [j]. Experimentally, the observable quantity is 
the mass difference (also named oscillation frequency) Am q , where the subscript q labels the light 
quark content (d or s) of the B meson. The mass difference is parametrized as [1] 



Am„ 



G 2 F m 2 w 



(1.1) 



where is the mass of the B 9 -meson, and V t * and V t t denote CKM matrix-elements. The Inami- 
Lim function, So [2], and the QCD coefficient, T]b UL can be computed perturbatively, while f q B Bq 
is the non-perturbative input: the decay constant f q and the B-meson bag parameter Bb ci - 

We define SU (3) breaking ratios as the ratio of a quantity for the B r meson over the same 
quantity for the Bj -meson. In particular we are interested in 



(1.2) 



and we also consider the ratio of the decay constants /bJ/b,,- Computing t, non-perturbatively 
allows one with additional experimental input to extract the ratio of CKM matrix elements 



Am s _ m Bs 
Am d 



m Bd \V td \ 2 ' 



(1.3) 



The phenomenological importance of B-B mixing is given by the fact that the ratio | V ; . s | 2 /| V t d\ 2 
constrains the apex of the unitarity triangle [3]. Experimentally, Amj and Am s are measured to bet- 
ter than a percent [4-7], whereas we know E, only to about 3%. Hence in order to get a stronger 
constraint on the apex of the unitarity triangle and consequently on new physics beyond the Stan- 
dard Model, we need to improve the determination of E, . 

The project presented here aims at this task by demonstrating the viability of our method 
on 16 3 x 32 x 16 dynamical domain- wall fermion ensembles generated by the RBC and UKQCD 
collaborations [8]. Our central values agree with the ones published in the literature, however 
our errors are large. A future follow-up project will hopefully lead to results with comparable 
uncertainties to those obtained by Gamiz et al. [9] and Evans et al. [10]. Details of the present 
computation will soon be published in [11]. 
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Figure 1. The dominant contribution to S-fi-mixing 
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# configs. 


ami 


ami, 


am x 


m n (MeV) 


APE HYP 


0.01 


0.04 


0.01, 0.0359 


400 


298 300 


0.02 


0.04 


0.02, 0.0359 


530 


298 300 


0.03 


0.04 


0.03, 0.0359 


630 


298 300 



Table 1. The light {mi) and heavy (/%) sea quark masses, the valence quark masses (m x ), the mass of the pions (m n ) 
and the number of configurations used in our simulations. 

2. Actions and parameters 

In order to accomplish the computation of fi-fi-mixing we choose domain-wall fermions [12, 
13] for the light quarks (u, d, s), the Iwasaki gauge action [14], and a static action with link- 
smearing for the ^-quarks [15]. 

Domain-wall fermions are formulated in five dimensions and posses an approximate chiral 
symmetry. Left-handed modes are bound to a 4-d brane at s = 0, while right-handed modes are 
bound to a 4-d brane at s = L s — 1. The overlap between both is exponentionally suppressed. 
Besides the approximate chiral symmetry, domain-wall fermions are advantageous because the 
renormalization is simplified due to reduced operator mixing. The choice of the Iwasaki gauge 
action is motivated by the fact that it improves chiral symmetry and reduces the residual quark 
mass in combination with domain-wall sea quarks [16] over e.g. the Wilson gauge action [17]. 

Our simulations use 16 3 x 32 x 16 dynamical 2+1 flavor domain-wall lattices [8] in which 
the up and down sea quarks are degenerate and the strange sea quark is slightly heavier than its 
physical value [18]. The domain-wall height is set to M=, = 1.8 and we have /3 = 2.13. We estimate 
the inverse lattice spacing to be aT x = 1.729(28) GeV and the residual quark mass to be am ves = 
0.00315 [18]. Hence a is approximately 0.11 fm and we have a (1.8 fm) 3 box. In Table |l| we list the 
values of the light (m/) and heavy (m,h) sea quark masses, the values of the valence quark masses 
(m x ), the mass of the pions, and the number of configurations used in our simulations. 

For the ^-quarks we use a static action which is based on the original lattice formulation by 
Eichten and Hill [15] but improved by link-smearing. This formulation corresponds to an effective 
action in the limit of infinitely heavy ^-quarks, i.e. corrections of order l//n/, are neglected. Ad- 
vantages of this static action are that the static quark propagator is given as product of gauge links, 
the propagator is 0(a) improved and has a simple continuum limit. We also enhance the signal- 
to-noise ratio by link-smearing in comparison to the original formulation [19]. For this work we 
utilize APE smearing [20, 21] and hypercubic blocking (HYP smearing) [22, 23]. 

In case of APE smearing with smearing parameter a = 1 , one computes all staples for a given 
link, adds them and projects the result finally back onto SU (3). HYP smearing consists of three 
steps of APE smearing restricted to the links within the hypercube surrounding a given link. For 
each step i we have a choice for the smearing parameters a,-. Following Delia Morte et al. [23] 
we choose (a,\, 012,053) = (1.0, 1.0, 0.5), which is commonly referred to as "HYP2". The required 
SU (3) projection is not unique: in the case of APE smearing it is performed by the unit circle 
projection method [24], while for HYP smearing we employ an iterative procedure which yields 
the SU (3) projection of a matrix V by seeking U m3X € SU (3) such that ReTr(?7 max \ / ^) is maximal 
[25]. Both projections are equivalent in the weak coupling limit [11, 24]. 
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smearing 


ca 




ZvA 


Z*sp 


APE 


0.0653 


0.9507 


0.7485 


-0.1448 


HYP 


0.1204 


0.9813 


0.8108 


-0.1448 



Table 2. Perturbative matching coefficients for bilinear and four-quark operators evaluated for APE- and HYP-smeared 
static-quark gauge links. 



3. Lattice calculation 

The perturbative renormalization of the heavy-light axial current and AB = 2 four fermion 
operator follow a two-step matching procedure. First we match the continuum QCD operators at a 
scale Hb in the MS-scheme using naive dimensional regularization onto operators in the continuum 
static effective theory at a scale ju. Next, we match the operators in the continuum static effective 
theory to operators on the lattice [26]. Combining the results of the two matching steps we obtain 
the perturbatively computed matching coefficients ca, Z®, Zva, and Z$p, which are listed in Table 
^| [27]. We account for the truncation of the perturbative series in our estimate of the systematic 
errors. 

For the bilinear operators we include the 0(a s pa) improvement term in our analysis 

Of" =Z*(1 +c A sinhK))^ at , (3.1) 

but we have not yet implemented O(pa) improvement for the four-quark operator. The renormal- 
ized decay amplitude <J>g n is related to the decay constant f Pq = ^™ I \J™%~ q m ^ * s obtained by 
computing on the lattice the ratio of two-point functions of the heavy-light axial current using local 
(L) and wall (W) sources/sinks: 

#-n»Jl . lvm{ ' M ■ 0.2) 

'»» » L ^ v ww (lilll)e -"i,t>-*.) 



where trig is the unphysical B^-meson rest mass. For the matrix elements containing the four- 
fermion operator we compute additional three-point functions ^ using box sources/sinks (B) [28] 
in case of APE smearing: 

M% = lim 2 ^ V/ ' ' ; (3.3) 

6 >f»t»<o ^ BB (t,tf)^ BB (t,t ) 

and wall sources/sinks (W) for HYP smearing: 

Because the four-quark operators of different chiralities mix under renormalization, we obtain the 
renormalized expression for the matrix element by 

= Z VA M l v % +AA + Z SP M^ +PP . (3.5) 
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. APE: <S> B J®B J = 1-214(79), } 2 /dof = 0.6 
l HYP: <6 B ,/9 Bd = 1.201(42), X 2 /Aal = 0.7 
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Figure 2. Chiral extrapolation of ^bJ^b,, = \J m B, l m B d ■ /b, I fB d (left) and ^/ms JiriBj ■ % (right). The APE data points 
are shown as blue triangles, whereas the HYP data points are shown as red squares. The color of the shaded (hatched) 
error bands match those of the APE (HYP) data points. The physical average u — d quark mass is indicated by the dashed 
line and the black dot denotes the physical strange quark mass. Only statistical errors are shown. 



Finally, we extrapolate our lattice data to the physical quark masses and the continuum us- 
ing next-to-leading order partially quenched SU (3) heavy-light meson chiral perturbation theory 
Schematically, the expressions for the SU (3) breaking ratios are given by 



1 + "chiral logs" + — — — yc va i(m s / — mi), (3.6) 



<D™ b (AnfY 



J = 1 + " cniral lo g s " + J^]f d ~™ 1 ( m *' " m/ ) ' (3 " 7) 

where the quark masses are expressed as dimensionless ratios and "chiral logs" denote non-analytic 
functions of the pseudo-Goldstone meson masses. Performing a linear, one-parameter fit of our data 
with respect to the expressions above we are able to extract the physical value for <J>J en /<J>^ en and <^ . 
As input parameters we use: A x = 1 GeV, = 2.35(16) [18], f n = 130.4 MeV [7], g B *Bn = 0.516 
[29], am U( i + am res = 0.001300(62), am s + am res = 0.0375(17) [18]. A detailed description of our 
chiral extrapolation will be included in [11], where we also will discuss the alternative of using 
heavy-light SU (2) chiral perturbation theory. 



4. Results and Conclusion 

We present the preliminary results of our chiral extrapolation in Fig. |[ The plot on the left 
shows the data and the fit for <t>B s /<i>B d , the plot on the right the outcome for ym]jJmB d • %■ In both 
cases only statistical errors are shown, which are computed following to Ref. [30]. In addition we 
estimate systematic errors by varying the input parameters around their uncertainty, considering 
a constrained linear fit as alternative fit function, and using power-counting for the discretization 
errors, the errors due to the renormalization factors, the finite volume errors I /nib corrections. All 
of these errors are listed in Table || and discussed in detail in [11]. 

Finally, we use the experimentally-measured ratio of the masses m B o/m B o = 5366.6/5279.5 = 

s d 

1.0165 [7] to obtain the following values for the SU (3) -breaking ratios of B-meson decay constants 



5 



Neutral B-meson mixing from dynamical lattice QCD 



Oliver Witzel 



hjh d % 



uncertainty 


APE 


t ivn 

HYr 


Arb 


T T \ ' T 1 

HYr 


statistics 


7% 


4% 


5% 


4% 


chiral extrapolation 


11% 


11% 


12% 


11% 


uncertainty in g B * B7 t 


4% 


4% 


3% 


3% 


discretization error 


3% 


3% 


4% 


4% 


renormalization factors 


0% 


0% 


2% 


2% 


scale and quark mass uncertainties 


1% 


1% 


1% 


1% 


finite volume error 


1% 


1% 


1% 


1% 


1 /ntb corrections 


2% 


2% 


2% 


2% 


total systematics 


12% 


12% 


13% 


12% 



Table 3. Total error budget for the SU (3)-breaking ratios fgj 'fg d and S, rounded to the nearest percentage. 



and mixing matrix elements: 



fa* 



1.20(08) (14) 
1.19(05)(14) 



APE 
HYP 



and 



1.19(06)(15) 
1.19(05)(14) 



APE 
HYP 



(4.1) 



where the first errors are statistical, the second are the sum of all systematic errors added in quadra- 
ture. Currently, we are updating our analysis and hence these values may change in our publication 
[11]. We find that both smearings used agree very well indicating that the discretization errors are 
small in the ratios. When comparing these new results to the ones in the literature published by the 
HPQCD collaboration [9] and presented by the FNAL-MILC collaboration at Lattice 2008 [10], 
we also find good agreement (see Fig. ||). However, our errors are large and we look forward to 
improve upon them in future works. 
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Figure 3. Comparison of lattice QCD results for / fg d (left) and t, (right) . Our new preliminary data points are marked 
as blue triangle (APE) and red square (HYP). The gray circle marks the values published by the HPQCD collaboration 
[9], the beige diamond the preliminary value presented by the FNAL-MILC collaboration at Lattice 2008 [10]. 
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